Spin-dependent processes at the crystalline Si-Si02 interface at high magnetic fields 
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An experimental study on the nature of spin-dependent excess charge carrier transitions at the interface be- 
tween (111) oriented phosphorous doped ([P]~ lO'^cm^-^) crystalline silicon and silicon dioxide at high mag- 
netic field {Bq ~ 8.5T) is presented. Electrically detected magnetic resonance (EDMR) spectra of the hyperfine 
split -"P donor electron transitions and paramagnetic interface defects were conducted at temperatures in the 
range 3 K< T <12 K. The results at these previously unattained (for EDMR) magnetic field strengths reveal 
the dominance of spin-dependent processes that differ from the previously well investigated recombination be- 
tween the donor and the ?(, state, which dominates at low magnetic fields. While magnetic resonant current 
responses due to ^ ' P and Pf, states are still present, they do not correlate and only the Pj, contribution can be asso- 
ciated with an interface process due to spin-dependent tunneling between energetically and physically adjacent 
Pb states. This work provides an experimental demonstration of spin-dependent tunneling between physically 
adjacent and identical electronic states as proposed by Kane for readout of donor qubits. 



PACS numbers: 76.30.-v, 71.55.Cn, 73.40.Qv 

Phosphorus doped crystalline silicon (c-Si:P) is one of 
the most widely utilized semiconductor materials, with ap- 
plications ranging from conventional microelectronics^ to 
proposed and presently widely investigated concepts for 
spintronics" and spin-based quantum information processing 
(QIP)1. Silicon based spin-QIP and spintronics concepts aim 
to utilize the comparatively weak spin-orbit coupling present 
in this material, and the correspondingly very long spin- 
coherence times^^, as well as the impact of spin-selection 
rules on electronic transitions which can be used for spin 
readout'*. Most of these applications involve electrical trans- 
port and spin manipulation at or near the silicon-silicon diox- 
ide (Si02) interface, making the understanding of spin pro- 
cesses in this region extremely important. Numerous studies 
of spin-dependent transport and recombination at the inter- 
face between c-Si:P and Si02 have recently been undertaken, 
with the aim of identifying and understanding these mecha- 
nism s^'^-^ , and showing that they can be utilized for the obser- 
vation of very small ensembles of donors** and coherent spin 
motion^'^. Additionally, spin dependent transport in two di- 
mensional electron gases at the c-Si/Si02 interface has been 
demonstrated"'-''. However, no systematic study of such pro- 
cesses at high magnetic field has been conducted to date, with 
the only data at magnetic fields Bo > 400 mT given by a single 
electrically detected magnetic resonance (EDMR) spectrum 
recorded at Bq = 7.1 T and a temperature T = 4 K 

In the following, a systematic investigation of the spin de- 
pendent processes at the interface between c-Si:P and Si02 
are presented for high magnetic fields (Bo ~ 8.5 T) at tem- 
peratures in the range 3 K< T <12 K.We show that the dom- 
inant spin-dependent recombination mechanism at low mag- 
netic fields, recombination between ^'P and Pt, centers, is not 



seen at high fields. Instead, transitions involving only ^'P 
donors or Pb centers dominate the observed EDMR signals. 
This study focuses in particular on the nature of the Pb only 
transition which has previously been observed at low mag- 
netic fields and nominally undoped c-Si-Si02 interfaces*', but 
has, however, not been observed in the presence of ^ ' P donors. 

Experimentally, we used prime grade Cz-grown c-Si(lll) 
with a phosphorus donor concentration [P]w lO'^cm^^. The 
sample was contacted by thermal evaporation of a 100 nm 
Al-layer after a surface clean and subsequent removal of the 
native oxide by wet treatment with hydrofluoric acid. Fol- 
lowing this procedure and the structuring of the sample con- 
tacts by a photolithographic lift-off procedure, a native Si02 
layer was formed on the surface between the contacts due to 
the exposure of the sample to air at room temperature. Simi- 
larly to previous studies of spin-dependent recombination and 
transport at low magnetic fields'^ '"* '^, we used EDMR to in- 
vestigate these processes. With this technique, the photocur- 
rent through a sample is monitored while electron spin res- 
onance is used to manipulate the spin of paramagnetic cen- 
ters involved in spin-dependent transitions. The latter are 
detected by measurement of currents which change from a 
constant offset value under spin resonance conditions'^. In 
order to perform spin resonance at Bq ~ 8.5 T, the quasi opti- 
cal 240 GHz heterodyne spectrometer facility of the National 
High Magnetic Field Laboratory in Tallahassee, Florida was 
used'^. A sample compatible to the geometric constraints of 
the Fabry-Perot (FP) resonator of the spectrometer was used 
for the experiments, and is shown schematically in figure[T] It 
consists of an approximately 330 ^m thick, 8 mm x 8 mm sil- 
icon substrate sandwiched between two 160 fim thick quartz 
slabs needed as antireflection coatings to allow the 240 GHz 
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radiation to be coupled into the silicon bulk. The electrical 
contacts to the device are a 100 jUm wide grid structure con- 
sisting of five 10 /xm wide interdigitated fingers, with 10 /xm 
separation between opposite fingers. The contact fingers were 
approximately 6 mm long, yet fingers belonging to the two 
opposite contacts overlapped by only 1mm. This geometry 
ensured that (i) the active region of the sample was located on 
the optical axis of the FP resonator such that the microwave 
field Bi was maximal and homogeneous throughout the active 
area; (ii) the external contacts of the sample, which were con- 
tacted with silver paste, were well outside the Bi field such 
that they could not distort the FP resonator modes; and (iii) 
due to the length of the contacts (almost 12 mm, stretching 
across the entire beam diameter), all metal structures within 
the beam diameter were aligned perpendicular to the polariza- 
tion of the Z?i field, reducing loss due to microwave absorp- 
tion. The high magnetic field, Bq, is aligned normal to the 
sample surface. The photocurrent needed for the EDMR ex- 
periments was induced by white light (cold light) generated 
by a xenon discharge lamp, filtered of its infrared component 
and coupled into the sample via an optical fiber. 

For the data acquisition, the microwave radiation was mod- 
ulated which allowed a lock-in detection of the magnetic res- 
onance induced current changes. The relative current change, 
AI /I, observed is shown in Fig. |2^ for T —3 K, 6 K and 12 
K. There are three resonances clearly visible which (as for 
all data presented in this study) were fit by Gaussian func- 
tions. The two resonances at the highest magnetic fields are 
separated by 4.2 mT, as expected from phosphorus donor elec- 
trons due to their hyperfine coupling to the donor nuclear spin. 
These resonances were used to calibrate the magnetic field 
axis (with the applied frequency set to 240 GHz) due to the 
drift in the superconducting magnet and the internal field due 
to the polarized electrons. The resonance at lower magnetic 
field, at a g-factor of ^ = 2.0014, is assigned to the Pj, inter- 
face defect (a silicon dangling bond)^ due to the agreement 
between the experimentally determined ^-factor and the ac- 
cepted literature value of g — 2.0014, for Bq parallel to the 
(111) direction. In contrast to experiments at lower magnetic 
fields, the Pb resonance here is well separated from the two 
phosphorus resonances, outside the field range that connects 
the two hyperfine peaks. This is expected as while the mag- 
netic field separation of the Pb and phosphorus resonances 
depends on the g-factor difference, the phosphorus hyper- 
fine splitting is constant (4.2 mT) for high magnetic fields 
(B>>4.2 mT). 

The data show that the peak intensity. A, (defined as the in- 
tegrated resonance lines = areas of the Gaussian fits) of the 
Pb resonance and the sum of the two hyperfine coupled ^^P 
resonances have no correlation. While at T = 12 K, the Pb 
resonance is « 6 times larger than the sum of the areas of 
the two ^^P resonances, it is smaller at T = 6 K and T = 3 
K. Moreover, at T = 6 K, the signs of the ^'P resonances are 
positive, in contrast to the sign of the Pb resonance, which 
is consistently negative at all measured temperatures. These 
observations are in stark contrast to low field EDMR, where 



the dominance of the ^'P-Pb pair mechanism causes a com- 
plete correlation between the intensities of Pb and the inten- 
sity sum of the two ^'P peaks^. Additional evidence that dif- 
ferent spin-dependent processes dominate at high magnetic 
fields is given by the current transient following a single, 
short microwave pulse. Fig. |2]3) shows such transients off- 
resonance, and on-resonance for the phosphorus and the Pb 
peaks. These transients were taken with a pulse length of 8 
/xs. A short microwave induced current which decays after 
approximately t — 30 fis is seen in all traces. Additionally, 
a decrease in the current (A/ < 0) is seen for both the phos- 
phorus and Pb resonances. The return of the current to the 
steady-state value can in both cases be fit with a simple ex- 
ponential decay. This is different to the more complex double 
exponential quenching/ enhancement expected for the ^'P— Pb 
pairs that are visible at low magnetic fields. Additionally, the 
time constants of the two exponential recoveries, Tp^ ~ 64(2) 
/xs and Tphos = 23(2) /is for the Pb and phosphorus respec- 
tively, are very different, further indicating that the two res- 
onances are not due to the same processes. It shall be noted 
that the data in Fig.|2b) shows that EDMR at Bq « 8.5 T leads 
to a significantly smaller ratio between microwave induced 
artifact currents and spin-dependent currents than that seen 
at low fields. This makes high field EDMR on silicon signifi- 
cantly more sensitive than X-band EDMR typically conducted 
at Bo ~ 340 mT. 

In addition to the spectra displayed in Fig.lH we have mea- 
sured EDMR for a number of other temperatures between 3 
K and 10 K. Figure [3] shows the value of the peak areas of 
each of the three resonance peaks, namely the Pb and both 
the high field and low field phosphorus. In order to obtain 
the maximum area each peak was fit after correction of the 
lock-in phase due to the different dynamic behaviors of the 
Pb and ^'P signals. The magnitudes of the two phosphorus 
resonances are very similar, as expected due to the negligible 
nuclear polarization. However, as forshadowed in Fig.|2] the 
form of the temperature dependence is unexpectedly different 
from low field EDMR experiments. From high to low temper- 
ature, the ^'P signals are initially negative and smaller than the 
Pb signal, becoming larger and positive between T w 8 K and 
T w 4.5 K, before again becoming negative, but with a larger 
magnitude, for temperatures below T « 4.5 K. We note that 
the sign and amplitude of the resonance at the lowest tempera- 
ture recorded agrees with the single spectra reported by Honig 
and Moroz'^. 

The experimental data presented allows us to exclude a 
number of mechanisms as the source for the observed signals 
and temperature dependencies. First, the signals observed are 
probably not bolometric effects due to resonant heating since 
this is expected to exhibit non-linear monotonically decreas- 
ing temperature dependencies which for both the ^'P and the 
Pb signals are not in agreement with the observed conductivity 
changes. It is possible that the low temperature (T < 6 K) ^'P 
signal has some bolometric component, however, Honig and 
Moroz'^ have assigned this to a spin-dependent neutral donor 
capture and reemission process, based on a spectra quanti- 



3 



tatively similar to that presented here. Hence, we conclude 
that the signals observed must be due to spin-dependent elec- 
tronic transport or recombination processes (except for ^'P 
in the range T < 6 K). Second, from the different magni- 
tudes, the different signs and the different temperature depen- 
dencies of the ^'P and the Pb signals we conclude that the 
3ip p^ interface recombination mechanism that dominates 
spin-dependent recombination rates at low magnetic fields is 
not responsible for the EDMR signals at high magnetic fields. 
Thus, the observed ^'P and the Pb resonances must be due to 
independent electronic processes. The ^'P enhancement sig- 
nal for 4.5 K< r < 8 K is not understood at this time in ab- 
sence of a theoretical framework describing this temperature 
behavior. This signal can not be attributed to bolometric ef- 
fects for the reasons stated above and because of its sign, as 
resonant sample heating is expected to cause a decrease of the 
conductivity. Thus, the source of the ^"P enhancement signal 
can not be attributed to an interface effect and consequently 
the only signal that is clearly due to an interface process is the 
Pb-only transition. 

We now consider the underlying process leading to the Pb 
resonance. As the Pb center is a paramagnetic deep interface 
state^, spin dependent electronic transitions are described by 
a two spin-1 /2 pair model^ -''-. Thus, the EDMR signal from 
the Pb-only transition can be due to: (i) Strongly coupled elec- 
tron pairs, such as the charged excited state P,^*, that decay 
spin-dependently into a charged ground state P,^ (this model 
was first suggested by Friedrich et al— for the low field Pb- 
only signal observed at the interface of intrinsic c-Si to Si02), 
(ii) tunneling or (iii) energy loss hopping between adjacent 
singly occupied Pb ground states with identical or different en- 
ergies, respectively. Fig. |3]shows that the temperature depen- 
dence of the Pb-only signal can be well fit with a linear func- 
tion without offset (at T = the EDMR signal Al = 0) up to 
r « 10 K. This behavior clearly contradicts (i) as the EDMR 
signal is not expected to vanish for small T in this model. 
Models (ii) and (iii) involve transitions between adjacent Pb- 
centers as illustrated for the tunneling case (model (ii)) in 
Fig. Ilh). Using Simmons-Taylor statistics'^(an extension of 
Shockley-Read statistics to an arbitrary defect distribution of 
states) it can be shown that at finite temperature and under illu- 
mination, the occupancy of the Pb defects is given by a Fermi 



distribution /(£■) = +exp — -^-j^^ j , about a quasi- 
Fermi energy Egf with ks the Boltzmann constant. The den- 
sity of filled Pb states close to the quasi-Fermi energy is thus 
the Pb density of states (DOS), S{E), times the fermi distri- 
bution, S{E)f{E) as plotted in Fig.B;)(i) for T = 10 K. Sim- 
ilarly, the density of unfilled Pj^ states is given by D{E){1 — 
f{E)) and plotted in Fig. |4};)(ii) for T = 10 K, where D{E) is 
the P^ DOS. The P^ DOS is identical to the Pb DOS except 
that it is offset along the energy axis by the positive correla- 
tion energy A associated with double occupancy of the defect^ 
(i.e., D{E) = S{E — A)) as is illustrated by the sketch of the 
DOS in the c-Si bandgap shown in Fig.|4j3). We note that for 
the low temperature range investigated here, the DOS is ef- 



fectively constant (i.e., S{Eqf) = Sqf and D{Eqf) = Dqf) as 
the thermal energy is small {kBT = 0.86 meV at 10 K) com- 
pared with the energy scale over which the DOS vary. Next, 
we assume that every Pb state interacts with P,^ states within 
some interaction radius, r, which is independent of tempera- 
ture. If we now consider only spin pairs whose energy levels 
are aligned (as expected for the tunneling model (ii)), we ob- 
tain a density Hp, given by 



Hp = k/SqfDqf J f{E){l -f{E))dE = ^T (1) 

with the constant ^ — nr^SgFDQFkB- From Eq.[Tl we see that 
the density of spin pairs is linear in T, with no pairs at T = K. 
As we anticipate a proportionality between the number of spin 
pairs and the EDMR signal, the energy-conserving model of 
tunneling between Pb pairs is in agreement with the observed 
temperature dependence of the the EDMR signals. When we 
consider energy loss hopping transitions (model (iii)), the spin 
pair density becomes 

np = nr'-J SgFfiE) Dqf{1 -f{E'))dE'dE - (2) 

in contrast to the experimental results. 

We note that the transition from a ^'P— Pb process^ at low 
fields to a Pb only mechanism at high fields may be explained 
by considering the underlying spin dynamics. The strength 
of an EDMR signal becomes weaker as the ratio ^ = ^ of 
the difference of the Larmor frequencies in a pair, Ao), to 
jBi drops below 1, with y being the gyromagnetic ratio' 
Hence, at low fields EDMR signals are dominated by maxi- 
mized ^'P— Pb signals as ^ > 1 while ^ ^ 1 for the Pb pairs 
with A(B <C 76 1 due to the almost identical Lande-factors of 
two Pb centers. At high fields, EDMR signals are dominated 
by the Pb transitions since ^ > 1 for both ^'P— Pb and Pb — Pb 
pairs while at the same time the Pb — Pb pair density is signif- 
icantly higher than the ^'P— Pb pair density. 

In conclusion, we have shown that EDMR on Si:P at the 
highest magnetic fields reported to date allows us to observe 
the influence of Pb centers and ^'P donor atoms on spin- 
dependent photocurrents. In contrast to low magnetic field 
EDMR, there is no intensity or transient correlation between 
these signals and, in contrast to the Pb signal, we find no evi- 
dence that the ^'P signals are due to interface processes. The 
intensity of the Pb signal increases linearly with temperature, 
vanishing as T — > 0, which is shown to match the properties 
of charge carrier tunneling between adjacent Pb states. This 
effect is expected only at high magnetic fields and does not 
contradict the dominance of the well investigated ^'P— Pb in- 
terface recombination process at low magnetic field. Finally, 
we point out that the spin-dependent tunneling demonstrated 
in this paper is analogous to the mechanism proposed by Kane 
for readout of solid state donor qubits^. Whilst previous at- 
tempts to investigate this mechanism have relied on remote 
charge detection of the transfer of an electron between clusters 
of donors^'' and even two single donors^', this work demon- 
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strates a spin dependent electronic tunneling transitions be- 
tween localized defect sites. Note that the ^'P— Pb mechanism 
that dominates at low magnetic fields does not demonstrate 
this effect, as energy in not conserved. 
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FIG. 1: A sketch of the sample used in these experiments. The sam- 
ple is fabricated on Si(l 1 1) doped with lO' 5 phosphorus donors/cm-', 
contacts are made to the sample using aluminum contacts, and the 
surface is covered with a native SiOa layer. A simple measurement 
circuit is also shown. The sketch is not to scale, although the indi- 
cated dimensions are accurate. 
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FIG. 2: a) Plots of the relative photocurrent changes as a function 
of the applied magnetic field, So at temperatures T = 3K, 6K and 
12K. The data was taken with a photocurrent / = 600 nA. The sig- 
nal was obtained by microwave chopping (500/xs pulse length, 1 
kHz shot repetition rate) and lock-in detection of the photocurrent. 
b) Plots of the photocurrent changes as a function of time follow- 
ing a microwave pulse with Tp = 8/is length for magnetic fields 
off-resonance, on-resonance with the low field ^'P peak, and on- 
resonance with the Pj, peak, at T = 12 K. The data was collected 
from the average of many transients measured with a 1 kHz shot rep- 
etition rate. 




FIG. 3: Plot of the integrated area A under each of the three reso- 
nances as a function of the temperature. The linear fit to the Pj, data 
between 3 K and 10 K shows excellent agreement. The data at 12 K 
was taken at a different illumination intensity. The line joining the 
data for the -"P resonances is a guide to the eye. 
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FIG. 4: a) Cartoon of spin-dependent tunneling between two adja- 
cent uncharged Pf, centers where the ground state energy of one Pb 
matches the Pj^ /Pj, charging energy of the other, (i) Initially, the 
pair is in a triplet state and tunneling is not possible, (ii) Once the 
spin of one Pj, state is changed the pair has singlet content and tun- 
neling (iii) is allowed, (iv) An excess charge carrier pair recombines 
as it discharges the two charged Pb states, b) Sketch of the DOS of 
Pb centers within the c-Si bandgap^ . Note that the DOS of P^ and 
P^ is assumed to be identical but shifted by the correlation energy 
A. c) Plot, for r = 10 K, of the number of (i) filled Pb states, (ii) 
filled Pjj^ states and (iii) tunneling pairs as a function of energy about 
the quasi-Fermi level Eqf. The sharper pair distribution for T = 5 
K is indicated by the solid line. The plot illustrates that there are 
fewer pairs of adjacent Pb and Pj^ with matching energies when T 
decreases. 



